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BACKGROUND OF THE INVENTION 

Field of the Invention 

[0003] The preferred embodiments of the present invention relate generally to acoustic well 
logging. More particularly, the preferred embodiments relate generally to determining the acoustic 
velocity (slowness) and frequency dispersion of acoustic waves in an earth formation. 

Background of the Invention 

[0004] In the art of acoustic logging, the formation properties of interest are one or more of the 
compressional wave velocity, shear wave velocity and Stanley wave velocity. These acoustic 
velocities are indicative of the formation's ability to hold and produce hydrocarbons. Operation of 
a typical acoustic tool may involve placing the tool in a well bore and firing one or more acoustic 
transmitters periodically, thus sending acoustic energy into the formation. The acoustic energy 
propagates along the formation wall in one or more propagation modes, e.g. compressional or 
shear wave modes. Receivers on the tool, spaced apart from the one or more transmitters and from 
each other, receive acoustic energy as the various waves move along the formation wall past the 
receivers. The amplitudes and arrival times of the various acoustic signals at the receivers are 
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indicative of the velocity of the acoustic waves (or slowness, being the inverse of the velocity) 
within the formation. 

[0005] Determining the acoustic velocity with early wireline acoustic logging tools involves an 
adaptation of data processing techniques used in seismic work. In particular, a method called 
semblance or coherence is used. United States Patent No. 4,594,691 to Kimball, et al, (hereinafter 
the '691 patent) is exemplary of related art wireline sonic acoustic logging tools that use this 
semblance/coherence measure for determining acoustic velocities within the formation. As 
exemplified in the '691 patent, determining the acoustic velocity using a coherence calculation is 
merely a determination of the extent two or more received waveforms resemble one another. The 
semblance/coherence determination itself, however, is not at all concerned with the actual 
formation properties; rather, the power of the semblance/coherence measure is running the 
calculation on portions of each received waveform that should correspond based on an estimated 
slowness. 

[0006] The semblance/coherence measure of the related art is run multiple times at multiple 
slowness values, and the slowness values where the waveforms show the best 
semblance/coherence are assumed to be the correct slowness values for the formation. However, 
some earth formations exhibit a property where the slowness of an acoustic wave is a function of 
its frequency. For example, a low frequency acoustic wave may have a slower acoustic velocity 
than a high frequency acoustic wave within the same formation. In some respects, anisotropic 
earth formations may exhibit this property. Moreover, some earth formations may actually skew 
the frequency of the acoustic signals as they propagate through the formation in addition to having 
different acoustic velocities for the various frequencies. Running a semblance/coherence measure 
to determine acoustic velocity in formations where frequency dispersion is present tends to mask 
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the actual formation slowness in the plot because of the dispersion characteristics. Figure 4 shows 
a related art time versus semblance plot from which the formation acoustic velocity may be, to 
some extent, determined. It is noted that in the related art time versus slowness plot of Figure 4 
(with the semblance value shown in isometric lines), it is not possible to determine the acoustic 
velocity as a function of frequency. 

[0007] Thus, what is needed in the art is a way to determine the acoustic velocity (slowness) of 
acoustic waves in formations that exhibit frequency dispersion characteristics and/or propagate 
acoustic waves at different speeds depending on the frequency. 

BRIEF SUMMARY OF SOME OF THE PREFERRED EMBODIMENTS 
[0008] The problems noted above are solved in large part by a signal processing technique for 
acoustic logging devices. In particular, a downhole tool, whether a wireline device, a logging- 
while-drilling device or a measuring-while-drilling device, has a transmitter and a plurality of 
receivers spaced apart from the transmitter and from each other. Acoustic energy is launched into 
the formation from the transmitter, and the receivers detect the acoustic energy as it propagates 
along the borehole wall and in the formation, the receivers creating a plurality of time domain 
signals. The time domain signals are preferably Fourier transformed to create a frequency domain 
representation of each received signal. 

[0009] Preferably, values from each frequency domain representation of the received signals at a 
selected frequency are used to create a correlation matrix. Eigenvectors and eigenvalues of the 
correlation matrix are determined with all the eigenvectors of the correlation matrix forming an 
orthogonal basis or space. Preferably, at least the one highest order eigenvector is removed to 
create a subspace, the at least one highest order eigenvector corresponding to the signal component 
in the overall received signals, as opposed to the noise. Thereafter, a plurality of test vectors are 
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applied to the subspace, with each test vector based on a different slowness estimation, to 
determine an objective function value indicative of the extent to which the test vector maps to or 
may be represented by the subspace. If the test vector maps to the subspace, the corresponding 
slowness value thus maps to the noise, indicating the estimated slowness is not correct. However, 
if the test vector does not significantly map to the subspace, this means that the test vector (and its 
estimated slowness) more closely represents the actual received signal The objective function 
calculated for each test vector is preferably plotted in a slowness versus frequency plot. The 
process is repeated for multiple slowness values within a single correlation matrix, and for multiple 
correlation matrices across the frequency spectrum of the received signals. What is preferably 
produced is a slowness versus frequency graph which shows the slowness of the formation as a 
function of frequency, and thus shows the dispersion of the acoustic velocities within the 
formation. 

[0010] The disclosed methods comprise a combination of features and advantages which enable 
them to overcome the deficiencies of the prior art devices. The various characteristics described 
above, as well as other features, will be readily apparent to those skilled in the art upon reading the 
following detailed description, and by referring to the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0011] For a detailed description of the preferred embodiments of the invention, reference will 
now be made to the accompanying drawings in which: 

[0012] Figure 1 shows an acoustic logging device constructed in accordance with the preferred 
embodiments; 

[0013] Figure 2A shows an exemplary set of received signals in an acoustic logging device; 
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[0014] Figure 2B shows the frequency domain representation of the received signals of Figure 2A; 
[0015] Figure 3 A shows an exemplary Cartesian coordinate space; 

[0016] Figure 3B shows the exemplary Cartesian coordinate space of Figure 3A with one principle 
component (Z axis) removed; 

[0017] Figure 3C shows the Cartesian coordinate space of Figure 3A with a vector completely 

defined by the missing principle component; 

[0018] Figure 4 shows a related art time versus semblance graph; and 

[0019] Figure 5 shows an exemplary plot of the objective function calculated in accordance with 
M the preferred embodiment. 

j! NOTATION AND NOMENCLATURE 

* 25 

'£] [0020] Certain terms are used throughout the following description and claims to refer to particular 
fl system components. This document does not intend to distinguish between components that differ 
ff j in name but not function. In the following discussion and in the claims, the terms "including" and 
Q "comprising" are used in an open-ended fashion, and thus should be interpreted to mean 
"including, but not limited to. . .". 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0021] Figure 1 shows an acoustic logging device 10 constructed in accordance with the preferred 
embodiments. In particular, Figure 1 shows a tool 10 disposed within a fluid filled borehole 12 at 
some distance below the surface 14. The tool 10 is preferably suspended within the borehole by a 
multi-conductor armored cable 16, and the cable 16 preferably couples to a surface computer (not 
shown). The tool 10 preferably comprises an acoustic transmitter 18 and a plurality of receivers 
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20A-D spaced apart from the transmitter, and also from each other. Although the preferred 
embodiment comprises four receivers, any number of receivers may be used without departing 
from the scope and spirit of the invention. 

[0022] Thus, the tool 1 0 of the preferred embodiment is a wireline device; however, the principles 
described herein may likewise be practiced in a measuring-while-drilling (MWD), logging-while- 
drilling (LWD) or any downhole measurement system. It must be understood that while the 
preferred embodiments of the present invention are described in the context of a wireline device, 
the description in this manner should in no way be construed as limiting the general structures and 
processing methods disclosed herein to wireline devices. 

[0023] Operation of the wireline tool 10 preferably comprises placing the tool within a borehole 
and allowing the tool to drop or otherwise move beyond the formation of interest. The tool 10 is 
preferably slowly raised within the borehole, and while being raised the transmitter 1 8 periodically 
fires inducing acoustic energy into the formation. As the acoustic energy propagates through the 
formation, each of the receivers 20A-D detect a portion of the acoustic energy and convert the 
acoustic energy into electrical signals. Determining the characteristics of the earth formation 
preferably involves analyzing the received signals to extract the information such as the acoustic 
velocity (or slowness, being the inverse of the velocity) generally, the acoustic velocity as a 
function of frequency of the signal, and the like. 

[0024] In the embodiment where tool 10 is a wireline device, surface computers (not shown) are 
responsible for extracting this information. In an embodiment where the transmitter 18 and 
receivers 20A-D are located on a LWD or MWD device, making these acoustic measurements 
while the drill bit cuts through earth formations, some or all of the data processing required to 
extract the desired information may be completed downhole. 
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[0025] Figure 2A shows an exemplary set of received signals for the four receivers of the tool 10 
of Figure 1 . The exemplary set of received signals have their amplitudes normalized, and thus no 
difference may be discerned in terms of amplitude from the closest received signal A to the more 
distant received signal D. However, it is seen that first movement of received signal D is later in 
time than the first movement of the received signal A for the closest to the transmitter. 
[0026] Figure 2B shows a corresponding set of frequency domain representations of the four 
received signals A-D of Figure 2 A, with, in this example, the majority of the frequency response 
between 20 and 30 kilo-Hertz (kHz). Before proceeding, it must be understood that the waveforms 
in Figure 2A and the corresponding frequency domain versions of Figure 2B are merely 
exemplary. Other received signals and corresponding frequency responses may be observed, 
depending on the frequency of the source as well as the acoustic properties of the earth formation. 
[0027] The signal processing technique of the preferred embodiment comprises calculating a series 
of correlation matrices using values or amplitudes in the frequency domain representations of the 
received signals along constant frequencies. A plurality of correlation matrices are made, with 
each correlation matrix based on a different frequency within the frequency domain representation 
of the received signals. Within each correlation matrix, eigenvectors or component functions and 
corresponding eigenvalues are determined, with the highest order eigenvectors and eigenvalues 
preferably corresponding to the received signals, and the lower order eigenvectors and eigenvalues 
corresponding to received noise. All the eigenvectors taken together form an orthogonal basis or 
space. Preferably, however, one or more of the higher order eigenvectors are removed with the 
remaining eigenvectors becoming an incomplete basis or subspace. Preferably, a series of test 
vectors are applied to each subspace to determine whether those vectors map to, or may be 
represented by, the subspace. If the test vector significantly maps to the subspace, then a slowness 
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embodied in the test vector corresponds to noise in the received signals. If, however, the test 
vector (and its embodied slowness) do not significantly map to the subspace, then the test vector 
maps to the signal component of the received signals. The following brief digression highlights the 
concept of mapping to determine whether the test vector may be represented by the subspace. 
[0028] Consider for purposes of explanation a vector V in three-dimensional Cartesian coordinate 
space as shown in Figure 3 A. The vector V points in a particular direction that has components, or 
projections, in each of the X, Y and Z coordinate axis. Thus, the vector V may be represented as 
follows: 

V=X x i+X y j+X z k (1) 
where X x , X y and X z are projections of the vector V intheJC Y and Z directions respectively, and 
where / , j and k are unit length vectors pointing in the X Y and Z directions respectively. The 
vector V may likewise be written as: 

V = \l x 7 <g> J + X y J ® j + X z k ® [f + j + k] (2) 

With the values X x 9 X y and X z being eigenvalues, the vectors i , j and k being the eigenvectors 
of the system, and where <8> is the tensor, outer, or dyadic product. 

[0029] The eigenvectors in this exemplary case represent unit length vectors that, when taken 
together, define an orthogonal basis or space. Stated otherwise, these eigenvectors may be used to 
completely define functions within the space. The eigenvectors of the preferred embodiments need 

not necessarily be straight lines, and preferably are an orthogonal set of vectors based on the 
corresponding points in the frequency domain version of the received signals. 
[0030] Assume that in the Cartesian coordinate system a particular eigenvector is removed from 
consideration, say eigenvector k . If the presence of all three eigenvectors was necessary to define 
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the overall space, a system short one or more components may be called a "subspace." If an 
attempt is made to represent the vector V in the newly defined subspace having only eigenvectors 
in the X and Y directions, the closest the subspace may come to representing the vector V is as 
follows: 



This partial representation is symbolically represented in Figure 3B. Thus, the vector V is not 
folly defined using only the subspace as the vector has some components (in the Z direction). 
Further, consider the special case where the vector V only has components in the Z direction as 
exemplified in Figure 3C. In this situation, the vector V does not map to the exemplary subspace. 
[0031] Referring again to Figure 2B, preferably a plurality of correlation matrices are created. For 
example, a correlation matrix is run for all the data points intersecting the 20 kHz axis. Likewise, a 
correlation matrix is calculated for all the points intersecting the 30 kHz axis. Thus, there is 
preferably a series of correlation matrices, with each individual correlation matrix calculated at a 
single frequency within the frequency response spectrum. For notational purposes, these 
correlation matrices will be referred to as C/. Each correlation matrix is then separable as follows: 



where Qj is a matrix containing the set of eigenvectors for the particular correlation matrix, Q T f is 

the transpose of the eigenvector matrix, and V y is a matrix containing the eigenvalues for the 
particular correlation matrix. The eigenvectors are: 



(3) 



C f =Q T f V f Q 



(4) 



(6) 



e 



n 
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where e is a particular eigenvector, and where n is the total number of eigenvalues. And for 
completeness, the eigenvalues are thus: 



K 


0 


0 


0 




0 


0 


0 


"-. 0 


0 


0 


K 



(5) 



where the X are the eigenvalues, and where n is the total number of eigenvalues. 
[0032] In the preferred embodiments, at least one of the higher order eigenvector is removed from 
consideration. With the eigenvector corresponding to the higher order eigenvalue removed, the 
remaining eigenvectors no longer define a complete space, but instead define a subspace. If the 
correlation matrix C/ was an M x M matrix, the subspace matrix, N/ } would be an (M-l) x M 
matrix if only one eigenvector is removed, would be an (M-2) x M matrix with two eigenvectors 
removed, and so on. Thus, the subspace eigenvectors are: 

: (7) 

■ 

e 

_ n— r _ 

where r is the number of eigenvectors removed. Removing the higher order eigenvectors 
corresponding to the signal of interest may be likened to Figures 3B, 3C that have the Z axis 
component of the three dimensional Cartesian space removed. 

[0033] In the preferred embodiments, the largest amplitude components of the received signals, 
and thus the largest eigenvalues and eigenvectors, correspond to acoustic signals of interest, such 
as the compressional wave, the shear wave or the Stanley wave. In creating the subspace N/ 9 
preferably the higher order eigenvectors, corresponding to signals of interest, are removed. 
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However, it need not necessarily be true that only higher order eigenvectors should be removed to 
create the subspace. It is possible that, depending on the acoustic tool used, other propagation 
modes could exhibit greater influence in the overall received signal than the signals of interest. For 
example, it is possible that tool mode waves may be predominant in a system, and in such a case, it 
would be preferred to not remove the eigenvector corresponding to the tool mode propagation. 
Likewise, it is possible that the casing mode wave could have a high order eigenvector in the 
overall correlation matrices, and preferably this eigenvector too would not be removed when 
creating the noise subspace. Further, though not preferred, it is possible to create the subspace to 
include the eigenvectors for the signals of interest, and to exclude the eigenvectors corresponding 
to noise or undesirable signals. In this case, eigenvectors associated with background noise, most 
likely the lower order eigenvectors, could be removed to create a subspace containing the 
eigenvectors of the desired signal components. Finally, it may be possible to pick and choose 
among the eigenvectors to remove them as desired to create the subspace, for example, removing 
higher order tool and casing mode waves, as well as lower order noise eigenvectors to create a 
subspace comprising only the eigenvectors associated with the desired signals. 
[0034] In the preferred embodiments, for each correlation matrix created a series of test vectors, 
each test vector given the notation W/ 9 are applied to the subspace matrix (the correlation matrix 
with removed eigenvectors) to determine the extent to which the subspace Nf may represent the test 
vector Wf. Preferably, the value of an objective function is calculated which is indicative of the 
extent to which a test vector may be represented by the subspace. In the preferred embodiments, 
when the test vector Wf maps to the subspace Nfi the objective function is small. When the test 
vector Wf does not map to the subspace Nfi and therefore does not map to the noise, the objective 
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function gets very large. Preferably, the objective function is calculated using the following 



equation: 



1 



(8) 



[0035] Each value calculated for the objective function is preferably plotted in a frequency versus 
slowness plot. In the preferred embodiments, the value of the objective function is shown in the 
graph in color with red and hot colors showing larger objective function values, and blue or cooler 
colors showing smaller objective function values. Figure 5 shows an exemplary plot of objective 
function values of the preferred embodiment with isometric lines showing areas where the 
objective function values are substantially the same. Notice from Figure 5 that in the frequency 
range spanning 20-30 kilohertz (which is the bulk of the frequency response indicated in the 
Frequency response graph) that the slowness is well defined and showing a slowness of 
approximately 175 jis/ft (region 100 in Figure 5). Notice also that in the frequency range of 
approximately 35 to approximately 45 kHz, that the slowness is approximately 125 (is/fl 
(region 102 in Figure 5). Finally, notice how in the frequencies around 60 kHz, that two slowness 
values are indicated (about 275 for region 104 and about 75 for region 106). Regions 104 and 106 
may exemplify an anisotropic characteristics of the formation at issue. Thus, it is seen that using 
the processing method of the preferred embodiments, it is possible to ascertain an acoustic velocity 
or slowness of an earth formation based on the frequency of the acoustic signal. Comparing the 
preferred graph of Figure 5 to the related art time versus slowness graph exemplified in Figure 4 
for a same set of gathered data, it is clear that in the time versus slowness graph of Figure 4 one 
cannot ascertain the slowness of acoustic signals as a function of frequencies, and indeed it is 
difficult to precisely determine any slowness value. 
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[0036] Before discussing the specifics of the test vector, it is helpful to understand that a 
propagation delay, represented by d ■ s (distance multiplied by the slowness) is a phase shift in the 
frequency domain. That is, 

f(tMf(co) (9) 

f(t-ds)^f(a>)e jds (10) 
where the double sided arrow in equations 9 and 10 is a Fourier transform. The Fourier transform 
of the function f(t) is simply its frequency domain version f( w) . However, the Fourier transform 

of f(t-ds) creates the frequency domain representation f( co) as well as a sinusoidal component e! ds . 

Considering a plurality of frequency domain received signals for four equally spaced receivers on a 
tool 10, the following represents the relationship between the Fourier transform of the received 



signals: 



f,(t)<r>f 2 (a>) (11) 



f 2 (t-ds)^f 2 (o>)e Jdi (12) 



f 3 (t-2ds)^f 3 (m)e' 2ds (13) 

f 4 (t-3ds)^f 4 (6))e m (14) 

where the distance between the transmitter and the first receiver (equation 1 1) is ignored, and the 
distance between each subsequent receiver is the same. Thus, the difference in the frequency 
response representation of the received signals is the phase relationship as a function of the 
distance between the receivers. The distance d between the receivers is known. 
[0037] Thus, the test vector preferably takes a form as follows: 
W f =[l e jds e j2ds e j3ds e^*] 
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where d is the distance between the receivers, s is the estimated slowness for the particular test 
vector, n is the total number of received signals and r is the number of removed eigenvectors. 
Thus, a plurality of test vectors are applied to the subspace N/ 9 each test vector preferably using a 
different estimated slowness of the formation. 

[0038] The above discussion is meant to be illustrative of the principles and various embodiments 
of the present invention. Numerous variations and modifications will become apparent to those 
skilled in the art once the above disclosure is fully appreciated. It is intended that the following 
claims be interpreted to embrace all such variations and modifications. 
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